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We study a four-level atomic scheme interacting with four lasers in a closed-loop configuration with 
a <0 (diamond) geometry. We investigate the influence of the laser phases on the steady state. We 
show that, depending on the phases and the decay characteristic, the system can exhibit a variety 
of behaviors, including population inversion and complete depletion of an atomic state. We explain 
the phenomena in terms of multi-photon interference. We compare our results with the phase- 
dependent phenomena in the double-A scheme, as studied in [Korsunsky and Kosachiov, Phys. Rev 
A 60, 4996 (1999)]. This investigation may be useful for developing non-linear optical devices, and 
for the spectroscopy and laser-cooling of alkali-earth atoms. 



INTRODUCTION 



Absorption and emission of monochromatic light in two-level atomic transitions are well-understood processes in 
quantum optics Their properties, however, can change drastically if transitions to a third atomic level have to 
be included. This is the case, for instance, in the A configuration where two stable states are coupled to a common 
excited state by laser fields, thus providing two excitation paths which can interfere. This interference lies at the 
heart of Coherent Population Trapping (CPT) Here, destructive interference between the transition amplitudes 
gives rise to a superposition of atomic states (dark state) that is decoupled from coherent radiation but populated by 
spontaneous emission. Consequently the atom becomes "trapped" in this coherent superposition. 

For configurations like the A scheme, the relative phase of the laser fields does not affect the steady-state dynamics, 
in the sense that there always exists a reference frame in which the Rabi frequencies are real. This is no longer 
fulfilled in closed-loop configurations i.e. when a set of atomic states is (quasi- ) resonantly coupled by laser 

fields, such that each state of the set is connected to any other via two different paths of coherent photon-scattering: 
In this case, the relative phase $ between the transitions determines the interference and hence critically influences 
the dynamics and steady state of the system |[ |^ . Previous studies of closed- loop configurations often featured 
double-A systems, where two stable or metastable states are -each- coupled to two common excited states |[ 0, ||- 
These works have shown a rich variety of non-linear optical phenomena. 

In this paper, we investigate the phase-dependent dynamics of a closed-loop configuration, consisting of four tran- 
sitions driven by lasers. One ground state is coupled in a V-type structure to two intermediate states, which are 
themselves coupled to a common excited state in a A-type structure. We label this system the (diamond) scheme. 
The steady-state of the scheme, like that of the double-A configuration, is a periodic function of the relative phase 
$ between the excitation paths, which contribute to the scattering between any initial and final state of the scheme. 
In particular, the steady state is determined by the concurrence of the phase-dependent Hamiltonian dynamics and 
the relaxation processes. We will show that, depending on $ and on the lifetimes of the intermediate states, the <) 
system can show a variety of behaviors including population inversion, CPT, and phase-dependent refractive indices. 
It is worth noting that the doublc-A and the <^ schemes are governed by the same Hamiltonian, but are characterized 
by different relaxation processes. This results in critical differences in the dynamics, which we will point out in our 
discussion. 

Excitation configurations like the <^ scheme have been investigated in the literature as a model for observing 
pressure- induced resonances and can be found for instance in experiments with gases of alkali-earth atoms which 
aim at optical frequency standards [llo| or at reaching the quantum degeneracy regime by all-optical means . Our 
investigation may contribute to the spectroscopy of these systems and to the development of new and efficient methods 
of laser-cooling. 

This paper is organized as follows. In section II we introduce the model and the basic equations. We discuss the 
properties of the system, and identify some relevant parameter regimes. In section HI we calculate the steady-state 
solutions as a function of the relative phase between excitation paths for certain parameters, and discuss the results. 
In Section IV we draw the conclusions and discuss some outlooks. 
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THE MODEL 

In this section we introduce the model, using a density-matrix formahsm. We move then to a reference frame where 
the phase dependence is exphcit, and formulate the optical Bloch equations. Finally, we discuss the symmetries of 
the system as a function of the phase and of the relaxation processes. 

The Master equation for a single atom 

We consider a gas of atoms of mass m. The atoms are free, and interact with a multi-chromatic light field. For 
a sufficiently-dilute gas, each atom interacts individually with the light, which couples to a set of atomic levels as 
depicted in Fig. 1. This set contains a ground state two intermediate states |2) and |3) and an excited state |4). The 
transitions |1) |2), |3) and |2), |3) — > |4) are optical dipoles with decay rates 72, 73, 742, and 743, respectively. Each 
dipole transition is driven resonantly by a laser, which is here considered to be a classical running wave, propagating 
along the z axis. The laser coupling to the transition \i) \j) is characterized by frequency ujij and wave vector 
kij, while the strength of the coupling is given by the Rabi frequency gye'-^^'J , where gij is real and Xij is a constant 
phase, determined by the phase of the atomic dipole and by the phase of the laser at time t = and position z — 0. 
The state of one atom at time t is described by the density matrix a, which obeys the Master equation: 

d X 

-g-^a=-[Hit),cT]+Ca. (1) 

Here, the Hamiltonian H{t) contains the coherent dynamics of the atom and the Liouvillian C describes the relaxation 
processes. The Hamiltonian depends explicitly on time and reads: 

Hit) = 

i=2,3 

-h<?j-4e"'^"^''*"'=^'*^+^^''^|4)01 4-H.c.). (2) 

The first term corresponds to the kinetic energy of the atomic center of mass, where pz is its momentum along the 
z axis (for simplicity, we consider only the component of the atomic motion along the z axis). The second term 
corresponds to the internal Hamiltonian, where hiOj denote the energies of the atomic states |j) relative to the energy 
of the state The remaining term describes the atom-laser interaction. 

The relaxation processes are assumed to be solely radiative, and the Liouvillian £ in (|^) has the form: 

Co ^ 74. /''^duAA4,(u)e-''=^^«-|j)(4|a|4)(j|e''=^^«^ - | (|4)(4|<7 + a|4)(4|) (3) 



J=2,3 

f-1 



E 



=2,3 



7-,- / du7Vji(u) 



-1 



where 74 = 742 + 743 is the total decay rate from level |4), and Nij{u) is the dipole pattern for the spontaneous 
emission of a photon on the transition \i) |j) (j = 2,3,4), i.e. the probability of emission at the angle 9 from 
the z axis, such that u = cos 9. This term describes the diffusion in the dynamics of the atomic motion due to the 
randomness of the incoherent events. In the limit in which the center-of-mass motion can be treated classically, and 
for sufficiently short times this effect can be neglected so that the explicit form of the dipole pattern is irrelevant. On 
this time scale and for a dilute gas, we can neglect relaxations due to collisions between the atoms. This is the regime 
that we are considering in the main body of the paper. 

Change of reference frame 

The Hamiltonian is explicitly time-dependent. For the configuration we consider here there exists in general 
no reference frame in which this explicit time-dependence can be eliminated. This is a characteristic of " closed- loop" 
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configurations Q, and it is a manifestation of the intrinsic phase sensitivity of the dynamics. However, in an adequate 
reference frame, the Rabi frequencies can be chosen such that only one is complex, with its phase $ being a function 
of all laser phases. Without loss of generality, we move to a reference frame where $ is associated with the laser 
coupling to the transition |3) — > |4), so that the coherent dynamics of the system is now described by the Hamiltonian 



n2 



4 



2m ^-^ 

+ ^ (ffi2 12> (1| + gi3|3> (1| + .92414) (2| + 5346^* |4) (3| + H.c.) , (4) 

where $ = z). The new Hamiltonian is connected with Hit) by the relation: ~ U^^ H{t)U — ifiU^^dU /dt, 

where U(t) is a unitary transformation, which reads 

U{t) = exp{-i[(tJi2 +C^24)<- (fcl2 +fc24)^+ (Xl2 +X24)] |4)(4|} (5) 

X exp {-i [{L0i2t - ki2Z + xi2)|2) (2| + (wiat - kirz + Xi3)|3)(3|]} . 
The detunings 5j in (^) are given by 



82 = {uj2 — ^12) + 
Sd. = (^^3 - W13) + 



TO 



2m ' 
2 

13 



TO 2m 

{ki2 + k24)Pz . Hki2 + fc24)^ 



S4 = {uji - W12 - W24) + . „ 

TO 2to 



and the phase $ is defined as 



where 



$ = Aujt - Akz + Ax, (6) 



Au = UJ12 + UJ24 - ^^13 - ^34, (7) 
Afc = ki2 + k2i - A:i3 - ^34, (8) 

Ax = Xl2 + X24 - Xl3 - X34- (9) 

The phase $ is the relative phase between the two excitation paths characterizing any transition between two atomic 
states. This phase is in general time- and position-dependent, and it results from the multi-photon detuning Aw, the 
wave-vector mismatch Afc, and from the initial laser and atomic-dipole phases. Ax. 

In the new reference frame the density matrix is p = U~^{t)aU{t), and its evolution is described by the Master 
Equation 

d 1 



(10) 

In the next section we derive the equations of motion for the matrix elements of p. 



Optical Bloch Equations 

In the following we assume a thermal distribution of atoms at temperature T, such that the thermal energy fc^T 
is much larger than the recoil energies h^kf^/2m. In this limit, we treat the atomic motion classically, and neglect 
the effects of the photon recoil on the center-of-mass dynamics . Then, the one-atom density matrix is given by 
P — J dzdpzn{z,pz)p{pz, z), where n{z,pz) is the atomic density as a function of the position and the momentum, 
and p{z,pz) obeys the Master equation at the parameters z,pz- Denoting with pij = {i\p{z,pz)\j) its elements, 
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with i,j = 1,2,3,4 (we drop the exphcit dependence on the parameters {z,pz) ), the Optical Bloch Equations (OBE) 
have the form: 



Pll = 

P22 = 

P33 = 

Paa = 

P12 = 



.512/ ^ , -913 , X 

1— (Pi2 - P21) + 1— (Pi3 - Pai) 

+I2P22 +73/O33, 

■ 912 f X , . 524 , X 

-1— (Pl2 - P2lj + 1— (P24 - P42j 

+742/344 - 72P22, 

.313/ X , . .934 / i$ _i$ N 

-1— (Pi3-P3i) + i— (e P34-e P43) 
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\5o 
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P24 
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.534 

1— e P24, 



(11) 

(12) 

(13) 
(14) 

(15) 
(16) 

(17) 

(18) 
(19) 

(20) 



and pji — (pij)* . In the following we will refer to the diagonal elements as "populations" , giving the occupation of the 
atomic states, and to the off-diagonal elements as "one-photon coherences" or "two-photon coherences" , depending 
on whether the involved states are coupled at lowest order by the scattering of one or two photons, respectively. 



Discussion 



Equations (|ll| - P0[) exhibit a parametric time dependence, which enters through the phase <I> as in (|^). Such behavior 
imposes limitations on the existence of a steady-state solution. Neglecting the coupling of the internal degrees of 
freedom with the external ones, two cases can be identified where the internal steady state exists: (i) when the scheme 
is driven well below saturation, and/or (ii) for Aw = 0. In (i), the processes leading to the absorption of two photons 
are negligible, and the relevant dynamics takes place between the ground state |1) and the intermediate states |2), |3), 
coupled in a (open-loop) V-type configuration |l^. In (ii), the steady-state solution is defined for any value of the 
other parameters. For Aw — and Afc = 0, $ is determined solely by the initial laser- and dipole-phases. While the 
dipole phases are fixed by the quantum numbers of the atomic transitions, the laser phases can be modified. In this 
regime, which we study in the remainder of the paper, the dynamics of the <^ configuration shares some analogies 
with an interferometer, the arms of which are formed by the multiphoton excitation paths and the phase difference 
between them is given by the relative phase Such analogy was first drawn in |^ for the coherent dynamics of a 
closed-loop scheme. In this spirit we interpret some of the results, presented below and in the next section, which 
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have been obtained in the presence of relaxation processes. 

An interesting manifestation of the phase-dependent dynamics is the probabihty of two-photon absorption on the 
transition |1) — > |4). In the hmit of weak excitations it has the form 



512524 , 513534 exp(i$) 



i72/2 



i73/2 



(21) 



Here, the first and second terms on the RHS describe the transitions via the intermediate states |2) and |3), respectively. 
From (^l]) it is evident that both the phase difference <I> as well as the ratio of the laser detunings 62, S3 and the decay 
rates 72,73 determine the interference between the two paths. In the special case of equal Rabi frequencies gij = g, 
and for S2 = S3, 'J2 — J3, the role of $ is singled out: 



Pi^4 cx cos" 



(22) 



Thus, the transition probability from the ground to the excited level is modulated by $. In particular, it is maximal 
for the values $ ~ 2mT (where n is an integer), while it vanishes for $ — {2n + 1)tt. At the latter value, no transition 
to |4) occurs. In the next section we will show that Pi^4 always vanishes for $ = (2n+ l)7r at steady state, even when 
the system is driven at saturation. We remark that the appearance of this behavior requires a "symmetric" excita- 
tion configuration, meaning that each two-photon excitation path from |1) to |4) has, separately, the same probability. 

A further understanding of the problem can be gained by moving to a suitable basis, following the analysis of [||. 
This basis is chosen appropriate to the structure of the Hamiltonian (^) and the relaxation processes in (||). For 
gij ~ g and Sj — the dynamics offers simple interpretations for $ = mr. 

We first focus on the values $ = {2n + 1)tt. Here, it is convenient to use the orthogonal basis set 
{|1),|4),|M/23(0)),|*23(^))}, where 



1 



with 9 = 0,7:. In this basis, the Hamiltonian (^) can be rewritten as 



H{{2n+l)TT) 



hg_ 
V2 



l*23(0))(l| + |4)(vl/23(7r)|+H.c. 



(23) 



(24) 



where we have omitted the atomic motion. Thus, ( p4[ ) describes two-level dynamics within the orthogonal subspaces 
{|1), |\['23(0))} and {|4), |vl'23(7r))}. These subspaces are coupled by spontaneous decay, and the coupling between 
states due to coherent and incoherent processes is represented in Fig. 2 a). From the structure of the decay, it is 
evident that the atom is eventually pumped into {|1), |\E'23(0))}. Hence, the steady state of the system for this value 
of the phase corresponds to that of the driven two-level transition |1) — > |4'23(0)). 

For $ = 2n7r we describe the system in the orthogonal basis set {|5'i4(0)), |\I'i4(7r)), [^'23(0)), |vE'23(7r))}, where 



1 

71 



ID 



with = 0,7r. In this basis, the Hamiltonian H can be written as 

H{2nn) = fig f|«'23(0))(*i4(0)| + H.c.l , 



(25) 



(26) 



and it describes a coherent two-level dynamics between the states |^'23(0)) and |5'i4(0)). The states |5'i4(7r)) and 
1^23(71")) are decoupled from the coherent drive because of destructive interference between the corresponding excita- 
tion paths, and from this point of view they are dark states. However, they are not stable, and decay with rates 74 
and 72 -I- 73 respectively. The level scheme in the new basis is plotted in Fig. 2 b). Here, it is evident that the system 
is incoherently pumped among the driven transition |^23(0)) |5'i4(0)) and the two dark states. One could say that 
the steady state is determined by the competition between the Hamiltonian dynamics and the relaxation processes. 
Thus, some localization in one of the dark superpositions (CPT) can occur, if this is more stable than the other, i.e. 
if the rate of pumping into it is much larger than its decay rate. In order to quantify this effect, we introduce the 
parameter a defined as: 



a 



74 

72 + 73 



(27) 
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Thus, for $ = and a ^ 1 (a <C 1) the dark state \^23{t^)) (|^i4(7''))) is long Uved with respect to all other states and, 
at steady state, it has a high probability of occupation. Such probability increases the more a differs from unity, and 
approaches 1 for a ^ oo (a ^ 0), corresponding to the system being trapped in |\I'23(7'')) (I*i4(7r)))- We will show 
that due to this effect, population inversion can occur on the transition |1) — > |2), |3) for a ^ 1 and on |2), |3) |4) 
for a <C 1. On the other hand, such behavior disappears as a approaches 1. For a = 1 and at saturation, the system 
is equally scattered among all states. 



Comparison with the double-A configuration 

In the absence of spontaneous decay, the <0> configuration is formally identical to the double-A scheme, extensively 
studied in the literature ^, Thus, the symmetries induced on the coherent dynamics by the phase are exactly the 
same We have discussed, however, that the steady state is critically determined by the concurrence between this 
symmetry and the relaxation processes. Thus, the introduction of the spontaneous decay leads to critical differences 
between the two systems. For an easier comparison we are labelling the atomic states of the double-A system as 
shown in the inset of Fig. 1. In this scheme the excited states |1) and |4) decay spontaneously into the stable or 
metastable states |2) and |3). In the scheme, the excited state |4) decays into the intermediate states |2) and 
1 3), which themselves decay into the ground state A first difference is that in the scheme the dynamics will 
be phase-sensitive only when the V- or the A scheme (or both) are driven at saturation, while below saturation it 
will reduce to the well-known V-configuration. In the double-A system, on the other hand, phase-sensitive dynamics 
survives also well below saturation 

When looking at the behavior as a function of $, the differences are more striking: at $ = {2n + l)7r, for instance, the 
(} scheme is pumped into the subspace |^23(0))}, which is a closed two-level transition, for the coherent drive as 
well as for the relaxation processes. In the double-A scheme, instead, the atom can be found in any of the four states 
due to incoherent coupling fl^ . 

At $ = 2mT, the role of the dark states differs between the two configurations. In the double-A system CPT occurs in 
the state 1^23(7'')), which is completely dark and, in the absence of other sources of decay, stable. In this configuration, 
and for equal decay rates from the excited states, the state |4'i4(7r)) is never accessed. In the scheme, instead, both 
dark states are accessed, and (partial) CPT occurs only when their decay rates differ substantially from one another. 



STEADY-STATE SOLUTIONS 

In this section we study the steady-state solution of (^^ as a function of the phase $. We consider laser frequencies 
and geometries that fulfill Ak — and Alu — 0, so that $ docs not depend on time and space. In order to obtain 
simple analytic solutions we consider resonant drives so that 6j = 0. Further, we assume that the moduli of the Rabi 
frequencies are all equal, gij = g, and that the decay rates fulfill the relation 72 = 73 = 7, 742 — 743 — 74/2. Under 
these assumptions, the system exhibits symmetry in the clockwise and anti-clockwise multiphoton excitation paths, 
the difference being the phase <I>. 

In this limit, we report and discuss the steady-state solutions of the OBE in ( pl|pO| ) as a function of the phase <I> 
and of the dimensionless parameters = 5/7 and a = 74/(27), defined in ([2^). We remark that in the following the 
rotated one-photon coherence 

P3i = P34e'* (28) 

is reported in the results. This frame allows to identify the real and imaginary parts of pi2, P13, /C24, Pai with the 
dispersive and absorptive response of the atomic medium to the fields which couple the corresponding transitions [ p^ . 

At the end of this section we will discuss experimental situations under which the assumptions given above are 
justified and discuss our results for generic values of 6, in relation to the assumption of classical motion, and for 
unequal Rabi frequencies and decay rates. 



Case a = 1 



We first consider the case a — 1. For convenience, we separate the real and imaginary part of the coherences, 
denoting them with ity = Ile{pij}, Vij = lm{pij}, respectively (here, U34 = Re{p34}, U34 = Im{p34}). The steady- 
state solutions of the OBE have the form: 
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where 
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:d 2 



2D 
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sin(f> 
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fl'' { Tsin^ i +27 
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(29) 

(30) 
(31) 
(32) 

(33) 

(34) 
(35) 
(36) 

(37) 

(38) 
(39) 



(40) 



The form of the solutions aUows to identify the contributions of the various multi-photon processes to the steady state. 
For instance, at second order in O (i.e. at second order in g/j) only pi4 depends on the phase while the populations, 
one-photon coherences and P23 are independent of <I>, and P44, P24, /534, w,i2! "la, and 7723 vanish. In fact, this limit 
corresponds to weak drives, and the relevant processes consist in resonant scattering on the transitions |1) |2), |3). 
Thus, at second order in O the system is equivalent to a V configuration driven below saturation. 

At higher order, the steady-state solutions arc phase dependent. This is evident, e.g., in the excited-state population, 
which is proportional to cos2($/2). In particular, at lowest order in fi, ^44 « ^l'^ cos2($/2). However, as ft is increased 
this modulated dependence of the populations is lost: at leading order in O, and for $ ^ (2n + l)7r, all states are 
equally populated. An exceptional behavior occurs at $ = {2n + 1)tt. Here, P44 = at all orders, while at leading 
order pn = 2p22 = 2/933 = 1/2- In Fig. 3 b) the populations are plotted as a function of the phase for = 2 and 
a = 1. (For comparison, Figs. 3 a) and 3 c) plot the populations for a <C 1 and a 1, respectively; we will discuss 
these regimes in the next subsection.) For the chosen parameters, pn, P44 vary with $, while P22, P33 are almost 
independent of the phase. 
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Fig. 3 e) and 3 h) show the one-photon coherences (|32-p5| ) as a function of $ for = 2 and a = 1. Their real parts 
vanish for <I> = titt, and one can easily verify from (|32|-35|) that this occurs at all orders in Q,. This is a feature of 
resonantly-driven two-level systems, and this result is consistent with the analysis of the previous section. Moreover, 
for $ — {2n + 1)tt one finds p24 — Pm — 0, which is consistent with the vanishing of P44. It is worth noting that U12 
and M13 can show additional zeros, as can be seen from their analytic form. These zeros depend on the values of 
and J7, which for a = 1 satisfy the relation sin2($/2) = 3(3 — 2^1'^) / SV,'^ . Thus, they exist only for a certain range of 
values of the Rabi frequency. Their existence can be interpreted as interference of multi-photon scattering at all orders. 

The two-photon coherence pu is proportional to cos($/2), a result in agreement with the calculation in (p^). The 
interpretation of the two-photon coherences becomes more transparent by employing the basis of the previous section. 
For instance, pu can be expressed as: 

Uli = i [(^'I4(0)|p|^'i4(0)) - (*i4W|p|*14(7r))] 

= Im{(*i4(7r)|p|«'i4(0))}. 

Analogue equations hold for U23 and V23. Thus, U14 = —1/2 (+1/2) corresponds to the system being in the state 
1^14(71')) (1^14(0))). The imaginary part V14 measures the coherence between these two states. We now look at 
( p6| - |39| ) as a function of <&, which are plotted in Fig. 3 k) and 3 n) for fl = 2. The behavior we observe is consistent 
with the above interpretation in the superposition basis, and with the discussion of the populations and one-photon 
coherences. At $ = 2mr the imaginary part U14 (^23) vanishes, supporting the hypothesis of no coherence between 
|5'i4(7r)) and |\E'i4(0)) (1*23(7'')) and |4'23(0))). Moreover, U14 = — M23 < 0, which imphes, after a straightforward 
calculation, that the probability to find the system in the dark states is 1/2. Thus, it is not proper to speak of CPT 
for these parameters. 

At leading order in Q the coherences vanish for $ — 2n7r, in agreement with the expectation that at saturation 
the system is equally distributed between all states. At $ = {2n + 1)tt one finds pi4 = V23 = while M23 is posi- 
tive and exhibits a local maximum. This is consistent with the picture of two- level dynamics between |1) and 1^*23(0)). 

So far we have discussed the case a = 1, when the relaxation rates of the two-photon coherences are the same. 
We have seen that the features of the phase-induced dynamics are always recognizable in the coherences. However, 
at steady state the atom is not localized in a particular atomic level or coherent superposition of atomic levels. In 
general, the dependence of the populations on the phase is washed out for increasing Rabi frequencies, except for the 
vanishing of P44 at $ = {2n+ l)7r. This is understood by considering that the steady state is given by the concurrence 
of the coherent drive, which has a phase-dependent symmetry, and the relaxation processes with a fixed structure of 
the coupling between the atomic states. For any value of $ 7^ (2?! + 1)tt, the two effects compete, and at saturation a 
kind of 'ergodicity' is recovered, so that the atomic states are equally populated. On the contrary, for $ = {2n + l)7r, 
an eigenspace of the coherent scattering processes exists and is preserved by the action of the incoherent processes. 
Consequently, at any value of H. and a the occupation of the state |4) vanishes. 

Case a 7^ 1 

Figures 3 a)-o) plot the steady-state solutions of the OBE for a = 0.1, 1, 10. Comparing the curves, we see some 
general features in the behavior at different a. For instance, the population of the state |4) is always zero for 
$ = (2n + 1)tt. This value of the phase is also a pole of the coherences p24, P34, Pi4, and of the real parts U12, U23- 
Here, the population of the state |1) and the real part of the two-photon coherence U23 exhibit a local maximum. 
These results are all consistent with the picture of two-level dynamics between |1) and 1^*23(0)), as discussed in the 
previous section. The steady-state values have a very transparent form for $ = {2n + 1)tt, and read: 

Pll = (41) 



^13 



1 


+ 2^22 


1 








1 




0, 









P22 = P33 = 1 , ,02 ' (42) 

(43) 
(44) 



1 4^2 ' 
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<2 = <3 = P24 = Pll = t^23 = Pll = 0, (46) 

which have been evaluated for O = 5/7 and an arbitrary value of a = 74/27. In these solutions, the parameter 74 
docs not appear, showing once again that the level |4) does not affect the steady-state dynamics for this value of the 
phase. 

A striking difference among the three regimes appears at values of the phase close to $ = 2n7r. Here, we find 

population inversion on the transitions |1) |2), |3) for a = 10, (|2), |3) |4) for a = 0.1), while the real part of the 
two-photon coherence U23 (^14) approaches the value —1/2. At this value of $, we write the steady-state solutions as 
a function of 7, a = 74/27 and Q = g/j: 

Pll = -^[a^(l + 2a) + fi2a(3 + 5a + 4a2) + fi4(i + 2a)], (47) 

P22 = P33 = ^[a{l + 2a) + n^a + 2)], (48) 

Pll = ^(l + 2a), (49) 

v!2 =vr3 = ^ [a(l + 2a) + 0^(2 + a)] (50) 

v'2l = v^l = ^(l + 2a) (51) 

<3 = ^[ail + 2a) + n\l-a)] (52) 

— [a{l + 2a) + n\l-a)] (53) 



tl4 



D 



*12 — "13 ~ '^14 — ^^23 ~ "24 — "34 



(54) 



where 

D = a^{l + 2a) + n^a{3 + 7a + 8a^) + 2n*{l + 4a + a^). (55) 
These results are plotted as a function of a in Fig. 4, by keeping 7 and as fixed parameters. 

Here, we see that for a <C 1 the system is localized in the atomic states |1), |4), and the coherence M14 has a maximum 
absolute value. In particular, for a ^ the populations of the states |2). |3) vanish together with the imaginary 
part of all coherences. In this case the atom is in the dark state |\l/i4(7r)), which is stable, and CPT occurs. Such 
localization persists for small values of a, although the populations of the intermediate states - and the incoherent 
scattering processes increase as a approaches 1 . It is interesting that for these (small) values of a the system exhibits 
population inversion on the transitions |2), |3) |4). Analogously, it can be verified that, fixed 74 and g, for 7 — > 
the system is trapped in the dark state |^'23(7r)) : CPT occurs in this coherence, and this implies population inversion 
on the transition |1) |2), |3). Note that the localization in an atomic superposition persists in the neighbourhood 
of the value of the phase $ = 2mT, as it is visible in Fig. 3 a) and c). For instance, for a 1 population inversion 
occurs on the transition |2),|3) — > |1) on an interval of values [2n7r — $Oj2n7r + $o]- The phase $0 satisfies the relation 
P22('^'o), P33('I'o) = Piii^o), and in general $0 can be said to separate two regimes, where the dynamics associated 
with the symmetry at phase 2mr or with $ = (2n + l)7r prevails. 

It is interesting to note that for $ = 2mr the populations and in particular the decay dependent population inversion 
show trends typical of a three- level cascade system while for $ = (2n + l)7r the system is effectively reduced to a 
V-configuration because of P44 = 0. 

Finally, we emphasize the additional poles of the one-photon coherences which we have identified in the analytical 

solutions for a = 1 . We have interpreted their origin as due to photon scattering at all orders. We remark that they 
appear in U12 and U13 for a = 0.1, see Fig. 3 d), and in U24 and U34 for a = 10, see Fig. 3 f). 



Discussion 



The analysis of this section is restricted to the choice of parameters Sj = 0, which corresponds here to values 
of Pz 0. This describes the behavior of a gas after Doppler cooling, at a thermal energy of kbT ~ ?ir/2, with 
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r = min(7,74). Provided that the hnewidth F is much larger than the recoil energies, T ^ h'^kf^/2m, the presented 
results describe sensibly the atomic response to the drive. 

When the medium is Doppler broadened, i.e. for kbT > fiV /2 (still keeping the constraint on the recoil energies), 
many features discussed for the case = survive and will appear in the signal measured over the ensemble, 
provided that 82 = 5^ = 5 5^ = Q, so that two-photon transitions are Doppler-free. This situation can be realised for 
degenerate intermediate-state energies (w2 = ^^3), resonant drives (tJi2 = 1^13 = 0^2 — t^i, tij24 = t^34 = 074 — 1^2), and 
laser geometries such that the wave vectors fulfill the relation fci2 — fcia ~ — fc24 — —^34 ~ k. In this way $ does 
not depend on time and space (Aw — 0, Afc — 0), and 5 is given by kpz/m. Also in this regime we find that for 
$ = {2n + 1)tt the population of |4) vanishes independently of 6, together with the coherences p24, P34 and P14. For 
$ ~ 2mr and a sufficiently different from unity, population inversion can be observed provided the atomic transitions 
are saturated [ p^ . 

Finally, we remark that only part of these considerations can be applicable to "asymmetric" configurations, i.e. for 
values of the Rabi frequencies, eigenenergies, decay rates, etc. , which change the structure of the Hamiltonian and 
relaxation processes, introducing thus either different weights to the interfering excitation paths, and/or additional 
relative phases, and/or different resonances. Here, the dependence of the steady state on the phase $ cannot often be 
simply singled out. The dynamics is a complex combination of all parameters, and exhibits an extremely rich variety 
of phenomena, which will be subject of future investigations. 

CONCLUSION 

We have studied the dynamics of a 4-level system interacting with lasers in a configuration which we have labelled 
the <0> scheme because of its geometry. This scheme has a closed-loop excitation structure |^ ^, i.e. any transition 
amplitude between two given states is the sum of two contributions, corresponding to two excitation paths, which may 
interfere. The dynamics is determined by a large number of parameters. Here, we have considered that both paths 
have the same weight, while they differ by a relative phase <&. We have discussed the origin of <i>, and investigated 
the steady state of the interacting system as a function of $, in the regime where the steady-state solution exists. 

For the chosen parameters, the steady-state solution is phase-sensitive. This is particularly evident in the coherences, 
whereas in general the phase dependence of the population is particularly enhanced for certain ranges of values of 
the relaxation rates. In particular, when the lifetimes of the intermediate states are considerably different from the 
one of the upper state, the system can exhibit population inversion for some values of the phase around <& — 2?i7r. 
We have interpreted and discussed this result in terms of coherent population trapping. Nevertheless, in all regimes 
here considered the population of the upper state vanishes for $ = {2n + l)7r. We have explained these behaviors 
using a convenient basis, showing that the dynamics is given by the concurrence of the Hamiltonian evolution, which 
is phase-sensitive, with the structure and non-unitarity of the relaxation processes. In particular, for $ = {2n + l)7r 
the steady state of the system corresponds to the steady state of a (closed) two-level transition. 

The phase dependence of the Hamiltonian evolution in closed-loop schemes shares many analogies with an atom 
interferometer |^ . The phase dependence survives also at steady state |^ , and the response of the system could be 
used as a device for measuring the relative phase between laser fields. For instance, in the <) system the phase could 
be measured through the population of the upper state. In fact, for sufficiently-weak fields, the functional behavior of 
this population is well approximated by cos^ ^/2, and the fluorescence signal from the upper state shows the features 
of an interference pattern which is sensitive to $. 

This study may be useful in the spectroscopy of alkali-earth atoms, currently investigated in experiments aiming at 
optical frequency standards |l^ or at Bose-Einstein condensation by all-optical means ||ll[. Further, efhcient laser- 
cooling schemes for these kind of atoms could be developed, by exploiting the phase properties due to the atomic 
motion in proper laser geometries [ pTf . 

Finally, the <0 scheme exemplifies a system where non-linear optics with resonant atoms can be realized. Here, the 
phase is a control parameter capable to change the response of the medium to the drive [|[ |[ 0, ||, |l8|, |l9| . This will 
be object of future investigations. 
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FIG. 1: Atomic level scheme and nomenclature or the relevant levels for the configuration. The inset displays the Double-A 
system for comparison. 




FIG. 2: Atomic level scheme in the basis of the states {|1), |4), |*23(0)), |*23(7r))} 
for the cases "I> = (2?! + 1)tt (a) and $ = 2nn (b). 
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FIG. 3: Steady-state values of populations and coherences as a function of the phase <I> for resonant drives and different decay 
rates: q = 0.1 with 74 = O.27 and g — 2^ (left column), a = 1 with 74 = 27 and g — 2^ (center column), a — 10 with 74 = 2O7 
and g = 74 (right column). Subplots a), b) and c): populations of pii (solid), p22 ~ paa (dotted) and p44 (dashed line). 
Subplots d), e) and f): real part of the one-photon coherences; U12 — — M13 (solid), and U24 — —U34 (dotted line). Subplots g), 
h) and i): imaginary part of the one-photon coherences; 1112 — V13 (solid), 1124 ~ V34 (dotted line). Subplots j), k) and 1): real 
part of the two- photon coherences; U23 (solid), U14 (dotted line). Subplot m), n) and o): imaginary part of the two-photon 
coherences; V23 (solid), «i4 (dotted line). Note the different scaling factor of the vertical axis in subplot o). 
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FIG. 4: Level populations and associated coherences for <1> — 2mi as a function of the balance of the decay rates a for g = 2^ 
(where 7 is kept constant), (a) populations: pn (solid line), p22 ~ P33 (dotted line), P44 (dashed line), (b) one-photon 
coherences: V12 = V13 (solid line), V24 = '534 (dotted line), (c) two-photon coherences: 1123 (solid line), 1114 (dotted line). All 
other coherences vanish. 



